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A simple system of two particles in a bidimensional configurational space S is 
studied. The possibility of breaking in S the time independent Schrodinger equation 
of the system into two separated one-dimensional one-body Schrodinger equations 
is assumed. In this paper, we focus on how the latter property is countered by 
imposing such boundary conditions as confinement in a limited region of S and/or 
restrictions on the joint coordinate probability density stemming from the sign- 
invariance condition of the relative coordinate (an impenetrability condition). Our 
investigation demonstrates the reducibility of the problem under scrutiny into that 
of a single particle living in a limited domain of its bidimensional configurational 
space. These general ideas are illustrated introducing the coordinates X c and x of 
the center of mass of two particles and of the associated relative motion, respectively. 
The effects of the confinement and the impenetrability are then analyzed by studying 
with the help of an appropriate Green’s function and the time evolution of the 
covariance of X c and x. Moreover, to calculate the state of the single particle 
constrained within a square, a rhombus, a triangle and a rectangle the Green’s 
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function expression in terms of Jacobi (^-function is applied. All the results are 
illustrated by examples. 


Keywords: Quantum boundary conditions, Confinement, Center of mass, Time evolu¬ 
tion, Jacobi ^-function. 

1. Introduction 

Over the decade, many papers appeared dedicated to the problem of quantum systems 
with boundaries. Several works are devoted to the problem of particles confined in a box, 
sometimes with moving walls m or specific shapes m ■ The confinement means the 
restriction on the motion of randomly moving particles, e.g. by the potential barriers. 

The confined systems have become a topical issue common to many research areas 
from condensed matter physics and quantum optics to biophysics. The first easy to un¬ 
derstand reason is that today an accurate quantitative prediction of the physical behavior 
of such systems is required by experimentalists since they cannot ignore confinement ef¬ 
fects when the spacial constraints stemming from the extreme miniaturization required 
for applications reach micro or nano-sizes. During the recent years some laboratory re¬ 
alization of new experimental systems like two and three dimensional graphite cones, 
carbon nano-tube rings [TOlITTj and torus shaped nano-rings iT2][l3] generated interest 
in the development of the idea of the quantum confinement. They can be used in devel¬ 
opment of nano and molecular electronic circuit devices. The second important reason 
is that the physical behavior of exemplary systems like a harmonic oscillator mns3 , an 
atom [IMS] or a small molecule, when subjected to such boundary conditions, may ex¬ 
hibit qualitative differences with respect to that available in the literature, for example, 
due to a breaking of the geometrical symmetry. This may lead to deep modifications on 
the eigensolutions of the system as well as to the need of a different way of treating the 
center of mass motion in the case of more than one particle in the system. 

In this paper, we first elucidate the breakdown of the separability of the center of 
mass motion from their relative motion for an unidimensional system of two, even non 
interacting, particles stemming from the confinement of the system. To this end, we 
first analyze the ” traditional” confinement constraint consisting on limiting the motion 
of the two particles inside a finite and the same interval / C 1R. In addition to such a 
boundary condition we distinguish the case when the relative coordinate is allowed to 
assume both positive and negative values from that when one of two particles is always 
on the same side with respect to the other one. We refer to this last situation speak¬ 
ing of an ’’impenetrability condition” describing it as an effective further constraint on 
the system. Recently unidimensional systems exhibiting this kind of restriction on the 
motion of the particles have been realized in laboratory j2'0] . Our target is to compare 
the separable motion of two totaly unconstrained particles, that is / = R and the lack 
of impenetrability, with the ones when at least one out of the two constraints is instead 
present. Our results clearly illustrate that the existence of a basis of factorized station¬ 
ary states of two even noninteracting quantum particles critically depends on whether 
and how the relevant dynamical variables get algebraically linked on the frontier of the 
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tridimensional domain out of which any wave function compatible with the assumed con¬ 
straints, vanish. In other words, we highlight that separability depends not only on the 
structure of the relative Schrodinger equation but also on the geometric shape of the 
normalization domain. This observation paves the way to the possibility of tracing back 
the confined motion of our unidimensional system of two noninteracting particles to that 
of a fictitious particle moving in a plane within a domain whose shape is determined 
by the restrictions imposed on our original two-particle system. The second part of this 
paper is thus dedicated to quantum billiards problems with plain domains as the square, 
the rhombus, the triangle and the rectangle. Also, the time evolution of the covariance of 
the center of mass coordinates is obtained. To this end, the theta-three Green’s function 
theory m is applied. It is worthy to mention that Jacobi 03 -function was used to study 
coherence states of a charge moving in constant magnetic field [22]. 

A quantum particle inside the potential of a special form is an interesting problem. 
There are many studies connected with the quantum and the classical properties of the 
bidimensional difficult geometries like, for example, the Robnik’s billiard and the Buni¬ 
movich’s stadium [231124] . the Sinai’s billiard [25]. The quantum particle in the triangle 
potential is studied in [26f[29| . The triangular shaped potential appears in different con¬ 
texts. For example, in the equipotential curves for the Henon-Heiles system i30]. Inside 
the boundary the potential behaves like the bidimensional harmonic oscillator, but in the 
billiard case, the particle has a free motion inside the triangular shaped domain. The 
brief review on the triangular billiard geometry problems is given, for example, in [31] . 
In this paper, we study four different shapes of the potential. It is shown how using the 
known boundary conditions for the particle confined in the square box with the side d , 
the boundary conditions for the boxes forming the rhombus, the triangle and the rect¬ 
angle are constructed. Moreover, the Green’s function using the Jacobi 03 -function for 
all four cases is obtained. 

The paper is organized as follows. In Sec. [A] the case of two noninteracting particles 
confined in the unidimensional box is considered. The boundary conditions are written 
for the center of mass motion. The case of the presence of the impenetrability condition 
is discussed. In Sec. m the problem of the single free particle motion on the plane in the 
domain of the motion being restricted by the unpenetrated walls forming the square, the 
rhombus, the triangle (the triangle billiard) and the rectangle are studied. In Sec. [4] the 
time evolution of the wave function is obtained. The time evolution of the covariance of 
the center of mass coordinates is shown for the bounded and the unbounded problems. In 
Sec. E] the Green’s function is used to find the time-dependent states for the covariance. 
The results obtained are illustrated on the example of the single particle motion confined 
in the square impenetrable box with the side d. 

2 . Two particles confined in a one-dimensional box 

Let us consider a system of two particles with masses nii and ni 2 confined in a 
one-dimensional box (tube) of a length d delimited by two impenetrable walls. The coor¬ 
dinates of the first and the second particle are 0 < x\ < d and 0 < X 2 < d, respectively. 
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The Schrodinger equation for the two noninteracting particles is the following 

h 2 / 1 d 2 1 d 2 \ 

-TUBf + =53) ,( * 1 '‘* )- “ ( * 1 '‘* ) - (1> 

It is known that in general the solution of m is the following 

^ free{xi,X 2 ) = {A sin(fcixi) + B cos{k\X\))(C sin(fc 2 x 2 ) + D cos(fc 2 x 2 )), (2) 

where fcj = \j ' lrr ^ E , i = 1,2. The four constants A,B,C,D must be determined taking 
into account normalization and the presence of the boundaries. In our case 4' f re e{x i, X 2 ) 
must vanish outside of the box and on the boundaries, that is 


W(xi = 0,x 2 ) = 4 f (xi,X2 = 0) = 0, ^(xi = d, x 2 ) = 4 '(xi,x 2 = d) = 0. (3) 

The solution 0 respecting such boundary conditions 0 can be hence rewritten as 


4' CO n/(xi,x 2 ) = A r sin(fciXi) sin(fc 2 ai 2 ), 


where ki = nni/d , fc 2 = 7rn 2 / d. N is the appropriate normalization constant. 
Introducing the well known center of mass transformation 


x = X\ — x 2 , X c 


TO1X1 + to 2 x 2 

M 


M = mi + m 2 , 


M = 


TOiTO 2 

M 


( 4 ) 


( 5 ) 


one can rewrite (JT|) as 


H 2 /I d 2 
2 \fi dx 2 


1 d 2 \ 

MdX 2 ) 


*{X c ,x) = E\b{X c ,x). 


( 6 ) 


For the problem of free particles this approach leads to the possibility of finding eigen- 
solution of (0 in the form of a product, e.g. 4'(X C , x) = \ki(X c )4' 2 (x). 

The representation of 0 in terms 0 is the following 


^confiAciX') 


X 


_ jy e ~iX c (k 1 +k 2 ) £ iX' 


n 2 k 1 —m 1 k 2 

- Tl - 


< e 2iX c k le -2ix^X _ 1 


)(< 


2iX c k 2 e 2ix^jp- _ 1 


where N is a normalization constant. 
From 0 we can write that 


( 7 ) 


Xl 


v , m 2 mi 

A c +—x, X2 = Xc -—x, 


0 < xi < d, 0 < x 2 < d. 


( 8 ) 


Hence, in the plane (X c , x) the domain of variability of X c and x turn out to be bounded 
by the four lines shown in Fig. [l] Moreover, we can deduce the geometrical boundary 
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Fig. 1: The domain of X c and x. 


conditions for the solution of © in the center of mass coordinates as 

$ (x c , x\X c + ^ x = d)= 0, T (X c , x\X c ~ r ^ x = d) = 0, (9) 

* (x c ,x\X c + ^x = O) = 0, * (x c ,x\X c - = 0) = 0. 

One should note that from the latter conditions the coordinates X c and x are algebraically 
related in the frontier of the domain. Hence, in the presence of the boundary conditions 
the separation of the variables is impossible [32] and we can not look for the solution of 
© as a factorized function of the two variables. Thus © can not be rewritten in two 
differential equations and relative boundary conditions each one depending only on one 
variable. We conclude that in the presence of the boundary conditions © the problem 
© is not separable. 

Let us search the solution of © in the form 

$(X x) = ^ e ^kl(aX c +bx) _|_ j^ e -ikl(aX c +bx)^^j e ik 2 (cX c +dx) _j_ jj e ~ik 2 (cX c +dx)^ 


Using the boundary conditions ©, we get o = c = 1, A = C = 1,B = D = — 1, 
b = — 7712 /M, d = mi/M and k\ = irni/d, k 2 = tt^/A Hence, the energy eigenvalues 
are 

h 2 kf H 2 kn 

En in2 = ^ + 7T^-- 

2777,1 2777-2 


Thus, the solution of © with respect to the new geometric boundary conditions © is 
the following 


^ni7i2 (X c , x') 


C( e ^i x c+^x) _ e -*=M* c +^z) )(e ^(* c -^) _ 



( 10 ) 








[Author and title] 


6 


that is 0 rewritten in the center of mass coordinates. The normalization constant can 
be found from normalization condition, rewritten in the notations (f5|) . Namely 


d 

/* 

/ 

/ dx 


J 

0 

J 

, mji 

\ rTf 


\]> 

r. 


l X) 



= 1 


and the constant may be assumed to be equal to A = If we assume that the crossection 
of the tubes possese a radius comparable with the liniar dimension of the particles, then 
they can not penetrate each other. Hence, we can think that the first particle is on the 
right hand side while the second one is on the left hand side of the box, e.g. x\ > x 2 - 
Using this additional condition and ® we can write 


.. m 2 „ m i 
X, + —x >X C - —x. 


( 11 ) 


Thus x > 0 that is always true for the present system, where 0 < x < d. That means 
that the additional condition on the penetrability of the particles influences the boundary 
conditions for the center of mass motion and the new domain is shown in FigJ^l Note 
that m does not fullfill the latter geometric boundary condition. In the presence of the 



Fig. 2: The domain of X c and x when the particles 
can not penetrate each other. 


inpenetrability condition, e.g. x\ > x 2 , but absence of the confinement (that means that 
x\,x 2 € (— 00 , 00 )) we can write (fill) and get the condition analogues to confinement 
from the one boundary. 

3. One particle confined in the two dimensional box 

The problem of the two particles confined in the one dimensional box is equvivalent 
to the problem of one particle confined in the two dimensional box. The Schrodinger 
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equation for such system is 

“i(^! + lg)® (s;i ' I2) = £ ' ,,(ll ' : ' 2) ' (12) 

Hence, in this section we study the the single free particle motion on the plane, but the 
domain of the motion being restricted by unpenetrated walls forming the the square with 
the side d , the rhombus with the side d, the triangle and the rectangle (see Fig. 0. 

Let us start from the square impenetrable box with the side d. The solution of the 
Schrodinger equation m is 


$( 21 , 2 : 2 ) = ^4sin(ni2i) sin(ri222), (13) 

where A is a normalization parameter. From the boundary conditions 
$(21 = ±d, 22) = $(21,22 = ±d) = 0 

we can estimate the parameters of the solution, ss. n\ = 1L ^ L . I n 2 = Ni,N 2 = 

0. : I- :2. 



Fig. 3: The domains of 21 , 22 . 


If we rotate the square box on 45° we will get the rhombus with the side d. As 
the solution of the Schrodinger equation m for such system we can take the following 
function 


$(21,22) = Bsin(m( 2 i+ 2 2 ))sin(n 2 ( 2 i — 22 )), (14) 

where B is a normalization parameter. The boundary conditions are 

$(21,22 = ±21 + y / 2 d ) = $(21,22 = ±21 — y / 2 d ) = 0 . 

Hence, the parameters of the solution (fT4l) are m = ^ 7 =^, 712 = N\, N 2 = 0, ±1, ±2,.... 

Moreover, dividing the rhombus domain at zero boundary, we get the triangle box. 
The latter case is called the quantum triangle billiard. As the initial step we study the 
billiard where the sides of the triangle are obtained by the intersections of the three lines 
described by the following equations 


21 = 0, 22 = 21 — V2 d. 


22 = —21 + V2 d 


(15) 
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Hence, the boundary conditions for such system are 

f(xi = 0 , 2 : 2 ) = f(xi,x 2 = ±xi + y/2d) = f{xi,X 2 = ±xi - y/2d) = 0. 

One should note here that all the latter conditions must be fulfilled simultaneously, which 
means that the solutions are equal to zero on the sides of the triangle given by (fl5l) . As 
the solution of the Schrodinger equation G 2 D that satisfy the latter conditions let us 
select the following function 


f(xi,x 2 ) = K($(xi,x 2 ) - $(- 2 : 1 , £ 2 )), 


(16) 


where K is a normalization parameter. The other parameters of the solution m are 

Finnaly, let us introduce the boundary conditions for the rectangle box. To this 
end we first study the square box with the side d and the Schrodinger equation of the 
following form 


~[hlk + ki g I *(«.«) = «■(«,* a ) 


(17) 


with the parameters of the Hamiltonian a and b. The energy is En i 

ni = n 2 = Ni,N 2 = 0,±1,±2,... and the solution of the latter equation 
is given by m■ Let us make the transformation to the new variables y\ = y/axi, 
y 2 = y/bx 2 - The transformed Schrodinger equation (1171) is the folowing 

1 / d 2 d 2 

2 \dyl + dy\ 

where the solution can be written as 


$ 


VL VL) = E( f> HL J'i 
Vb) y^Vb 


$ fiL 
Vv^’ Vb 


a ■ 1 y 1 \ ■ y 2 

= dsin ni~ 7 = sm n 2 —p 
yfaj V "Vb 


Thus, the boundary conditions for the rectangle domain with the sides d^/a 1 dy/b are 

(^yi = ±dVa 7 ^=\ = $ ^ Vl 


$ 


V 2 _ 

Vb 


,2/2 = ±dVb] = 0. 

y/a J 


We can conclude, that one can obtained the boundary conditions for the specific forms 
of boxes using the knowledge about the boundary conditions of the well known squar 
box case. 

4. Time evolution of the center-of-mass covariance 

Since the center of mass coordinates are dependent we are interested in there covari¬ 
ance time evolution for the bounded and unbounded problems. By definition, the square 
modulus of the wave function is the probability density function 


f(X c ,x,t) = $(X Ci x,t) 




[Author and title] 


9 


Thus, the expectations for X c , x and X c • x are the following 

oo oo oo oo 

E(a;) = f dx [ Xf(X,x,t)dX, E(X C ) = f dx f xf(X c ,x,t)dX, 


—oo —oo 


—oo —oo 


oo oo 


E(X C • x) = dx Xxf(X,x,t)dX 


—oo —oo 


and the covariance of X r and x is 


Cov(X c ,x,t ) = E(X C ■ x) - E(X c )E(x) = {$(X c ,x,t)\X c x\$(X c ,x,t)) (18) 

- ($(X C , t)\X C \$(X C , x , t)){$(X c , x , i) |x|$(X c , x, t)). 

However, we can just substitute the definition of X c and x in (1181) and get the equivalent 
result 


Cov(X c , x , t) 


E((mi£i + H 2 X 2 )(xi — x 2 )) — E(fj,ixi + M 2 £ 2 )E(aq — X 2 ) (19) 

Mi (E(a:f) - E 2 (aq)) - M 2 (E(x|) - E 2 (a; 2 )) 

Ml ({$(x 1 ,X2,t)\xl\${x 1 ,X 2 ,t)) - (($(x 1 ,X2:t)\x 1 \$(x 1 ,X 2 ,t))) 2 ^ 
M2 (['&{x 1 ,X2,t)\xl\$(x 1 ,X2,t)) - (($(x 1 ,X2,t)\x 2 \$(X c ,X,t))) 2 ^ . 


In the case of the equal masses the covariance is zero. Note that we must change the 
variables in the state function before integrating (fl9l) . 


4.1. Examples 

To obtain the covariance, one need to select some time state <k(X c , x, t ) ($(X C , x, t)). 
Using the evolution operator we can obtain any state |4>(f)) from the initial one |4»(0)} 
as 


| $(X c ,x,t)) = e |4>(X c ,ai, 0)). 
The evolution operator can be written as 


iHt 


_ L J—f — -L f 

ft O A/T L ft O ., 1 


g h = g h 2 M g h 2p, 


If the problem is unbounded the well known Baker-Hausdorff formula 
e ^Be~^ = B + £{AB} + £[A[AB]] + ... and 


can be used 


iHt 

X r xe 


= ( e S5W Ue h ™ 


’)(■ 


i£_f —±P—f 

e h2 » xe h2 ^ 


Xc + Ji t )( x + J'^ (20) 
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holds. Hence, the expectations are 

<$(i)|X c a:|S(i)> = (m\e^X c xe-^\m) = <*(0)| (x c + |$(0)), 

($(t)ix c |$ W ) = (m\e^x c e-^\m) = (m\{x c + ^\m), 

($(t)|x|$(f)) = ($( 0 )|e _?r xe __?r |$( 0 )) = ($( 0 )| + — t^j |$( 0 )) 

and the covariance of X c and x for the unbounded problem evolve in time as 

Cov(X c ,x,t) = ($(0)|^ e +-^tV®+^|$(0)) ( 21 ) 

- ($( 0)1 (x c + |$( 0 ))($( 0 )| |$( 0 )). 

For the bounded problem we can not use the Baker-Hausdorff formula since the 
variables X c and x are dependent. As an example, let us take as the state $(A C , x , t) the 
function that depends only from {ni,ri 2 } = {( 1 , 1 ), ( 2 , 2 )} and has the following view 

*{X c ,x,t) = A 1 ^ 11 {X c ,x)e- iE ^ h + A 2 ^ 22 (X c ,x)e~ iE ^ h . 

Hence, the probability density function is 

f(X c ,x,t) = A 2 * 2 1 (X c ,x) + A 2 * 2 2 (X c ,x) 

+2AiA 2 '&n(X c , x)^ 22 {X c , x) cos ((£n — E 22 )t/h ). 

From the normalization condition we can choose A\ = A 2 = §. Since for the bounded 
problem the coordinates X c and x are in the domain shown in Fig. [TJ the covariance 
m is the following 

Cov(X c , x, t) = - ^ 0 ^ 33 ^ 4 ^ ^(668288 - 61440tt 2 ) cos ((.Eh - E 22 )t/h) 

+ 102400 cos 2 {{E n - E 22 )t/h) + 8640tt 2 - 4608tt 4 + 50625 

o *r 2 2 . 

where En — E 22 = 2t p , /r = mim 2 /M. Certainly, the latter example is given for the 
illustrative purposes. To get the more general covariance we need to find the state |$(i)) 
for the bounded problem. To this end, the Green’s function theory can be applied. 

5. Green’s function and time evolution 

In [37j the propagator is investigated for the case of the one particle confined in the 
box. The relation between the wave function at time t and 0 is given by the following 
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formula 


oo 

'&(X c ,x,t)=ih J ty(xi,X 2 )G(X c ,x,xi,X 2 ,t)dxidx 2 , 


where G(X c ,x,Xi,X2,t) is a time dependent Green’s function or the propagator. In [35 
the Green’s function is represented in terms of theta-three function 

G(tvx) = ^( 03 (£, t ) -6 3 (v,t)), 
where the theta-three function is 

OO 

9 3 (C,r) = e 2inC e i7rn2r . 

n =—oo 

In the Whittaker and Watsou notation [ 36 j it can be represented as 

OO 

03 (C, 9) = l + 2 ^cos( 2 nC) 9 ” 2 , q = e™ T . (22) 

n= 1 

Hence, the propagator for the one particle is 

oo 

G(x,x',t) = — e 2in ( x - x ' ) e l ™ 2t . 


n ——oo 


It is obvious, that the propagator for the two particles can be written as 

G(X c ,x,x 1 ,x 2 ,t) = j-p (9 3 (X c ,t) - 9 3 (x 1 ,t))(6 3 (x,t) - 0 3 (x 2 ,t)) 


M 2 I E 


^ 2 inX c giTrnt 


E 


g2inxi ginn t 


\n— — oo 


n ——oo 


^2ikx^i7vk 2 t ^ ^ 


g 2ikx2 gi7tk z t 


2 ikx ink 2 t 

s e 

\k =—oo k= — oo / 

As the initial state ^(xi, X 2 , 0) one can select the following function 


$(xi,x 2 ,0) = ±- 


1 


f sin xjJ sin [ 


\/2d V ^dJ \d 




27 T 


sin —x i I sin — £2 


27 T 
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Thus, we can write 


d d 


9(X c ,x,t) = | 0 3 (X c ,t)O 3 (x,t) I I V(xi,x2,0)dxidx2 


0 0 


d d 


Oz(x,t) J J 'ff(xi,x 2 , 0 ) 0 3 (xi,t)dxidx 2 
0 0 

d d 

-0 3 (X c ,t) J J 9(xi,X2,0)03(X2 9 t)dXidX2 

0 0 

d d \ 

j j ^(xi,X2 : 0)6 3 (xi,t)0 3 (x2,t)dx 1 dx2j 


0 0 
ihy /2 
4 d 


2 °° 1 4- p 2idn 2 

e 3 (x c ,t)e 3 (x,t)—-(e 3 (x c ,t) + 6 3 (x,t)) £ + ~ i ™ t 


+ 


2 / oo 2 id n i i 

7r / e +1 


E 


_2 / °° ~2idn 


\n =—oo 


7T 2 — 4d 2 ?r 2 


E 


\n=—oo 


gZidn _ ^ 

7r 2 — dPri 2 


Using that 


we can write 


1 

ix 


/'•OO 

/ exp (—is#) dz 

Jo 

pOO 

/ exp (—iz(7r 2 — 4d 2 n 2 )) ck. 
Jo 


i{ 7T 2 — 4 d 2 n 2 ) J 0 

Substituting 'F(X c ,x,f) in (1181) the covariance time evolution can be obtained. 


5.1. The Green’s function and the special boundary shapes 

Let us find the Green’s function for the four boundary shapes presented in Sec. [3] 
By definition the Green’s function is the following 

OO OO 

G(xi,x 2 ,x' 1 ,x' 2 ,t) = ^2 ^2 ®N 1 ,N 2 {xi,X 2 ,t)<f>* Nl N2 (x[,x' 2 )exTp(-iE Nlt N 2 t). (23) 

Nl= — OO N2 = — 00 

For the square box domain the energy is En 1 ,n 2 = '' 1 2 ^ I n ' 2 = ^ 2 mE 2 ^ ’ ni = iL j 1 J 
7i2 = N\,N 2 = 0, ±1,±2,.... Hence, for the latter boundary shape the Green’s 
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function (l 23 l) is the following 

OO 

G(xi,X2,x' 1 ,x' 2 ,t) = A 2 sin(nixi) sin(n2X2) sinlniXi) sin(n2x' 2 )e~ tENlN2t 

N ±, N 2 =—oo 


A 2 ~ 

- 4 - E 


N i=—oo 


COS 


7 rfVi .A ( 7 riVi /,ii _j 

——(a-i + aq) - cos —— (aq - aq) e 


E 

N<2,= — OC 


COS 


TtN2 


nNo , , „ \ \ i7r2 N 2 , 


(X2 + X 2 ) ) — COS ( (x 2 — ^ 2 ) ) ) e 2md2 


where we used the known formula sin a sin (3 = ^ (cos(a — j 3 ) — cos(a + / 3 )). Using the 
following notations 

7T 7T 7T 

Vl = - *l), ^2 = 2^(351 + %), ?? 3 = —(x 2 - X' 2 ), 

7T , / 7Tt 

^ = 55 fe+ a:2 ) , T = ss; . 
we can rewrite the latter Green’s function in a short form 


A 2 


G(m,V2,V3,V4,r) = — (cos( 2 i ll N 1 ) - cos (2772-/^)) e v 


rNf 


Ni=—oo 


(cos(2? ?3 iV2) — 0x3(2774-^2)) e : 

N2 = — 00 


i-KTNi 


Since iVi = 0 , ± 1 , ± 2 ,... we can split the latter sums into two sums as 


A 2 


G(v 1 ,V2,V3,V4,t) = — J 2 (cos{ 2 r] 1 N 1 ) - 003(2^4)) e z 


rNf 


\N i = l,2,. 
— 00 


y (cos( 2 r] 1 N 1 ) - cos( 277 2 iVi)) e* 


rJVf 


^ {cos( 2 t] 3 N 2 ) - cos( 277 4 JV 2 )) e 

, iV 2 = 1,2,... 


inrN% 


y, (cos(2t7 3 -/V 2 ) - cos(2t7 4 -/V 2 )) e 1 


rJVJ 


JV 2 = -1,-2,. 


Using the definition of the theta three function (E^l) and the notation g = exp(77rr) we 
can finally write that the Green’s function for the square box with the side d is 
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For the rhombus box we can use the similar technics as in the previous case. The 
energy is E NliNa = , m = n 2 = JVi.iVa = 0,±1,±2,... 

and the Green’s function for the rhombus with the side d is the following 

G r ('ni,V2,V3,V‘i,q) = b 2 (6> 3 (Ci, 9 ) - 0 3 (C 2 ,q)) ( 03 ( 4 : 9 ) - 03 « 4 ,<?)), 
where we used the notations 

<i = urd (xi ~ x, ~ xt ~ x ^ 6 = 2 ^ <I,+I2+a: ' +l2) ' 

4 = ^=^1 - ^2 - 4 + 4 ), C 4 = - *2 + 4 - 4 ). 1 = e l7rT - 


For the triangle billiard the energy is equal to En 1: n 2 = 2 ” 1 ^' i ' 2 = - - , «i = 

v^d’ 712 = v^d’ Ni,N 2 = 0? =^=2,..., q = exp(2*7rr) and the corresponding Green’s 

function 

G b (rj 1 ,r]2,V3,Vi,Q) = K 2 ^ (0 3 (4: q) ~ 03 ( 4 : q)) (03 (C3,9) - 03 ( 4 , 9 )) 

- (03(-Ci,9)-03(-C2,g))(0 3 (-C3,9)-0 3 (-C4,«)) 


holds. 

Finally, it is easy to verify, that for the rectangle domain the Green’s function is 

G6(4:4:4:4:9a:9&) = A 2 (03(4) 9a) - 03(4:9a)) (03(4) 0b) ~ 03(4:9b)) : 

where the following notations were introduced 


4 = 

2^d lSl *'>’ 

4 = 

2^d ( " 1 + i '‘ ) - 

—irt 

Ta ~ 2 ad 2 1 

—nt 
n - 2 bd 2 

4 = 

2^d lB »>’ 

4 = 

■2V-ad (V2+V2h 

q a = e i7rr “, 

q b = e™ T \ 


5.1.1. Example 

As an example let us check the obtained Green’s function corresponding to the square 
impenetrable box with the side d. If the initial state is m one can write 


y(xi,x 2 ,t) = j J Asm(n 1 y 1 )sm(n 2 y 2 ) : ^ (o 3 (^(xi-yi),q) - 9 3 (^(x 1 + yi),q)j 

—d —d 

■ (^3(^(x2-y 2 ),q)~0 3 (^{x2 + y2),q)^dy 1 dy 2 (24) 
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= 4- / sin (™Wi) Q N h cos (-j-( x i-yi))-cos(^^-(x 1 + y 1 )))dy 1 

-d v >-' 

d oo 

• J sin (7123/2) q N H cos (-j-( x 2 - 2/2)) - cos(^^(x2 + V2)))dy 2 


N 2 = l 

The integrals are the following 
d 


f . /7rlVl \ ,-/Vl7r \ \ /■^ r i 7r / , 

/ Sm 1 — yi ) -y 1)) -cos(—— (xi +2/i)))dyi 

—d 


7riVi 


(27riVi — sin(27rTVi)) sin ( -p-x 1 ) . 


Hence, substituting the latter result in (1M1) we can write 

' ttN\ 


X fr(xi,x 2l t) = A 3 d 2 5^ 

N 1 , N 2 = 1 

Since A = 1/d we get the known state 


sm 


£1 sm 


nN 2 


-a:2 ■ 


^(xi,x 2 ,^) = a 

Ai,A 2 = 1 


( iV 2 +iv 2) 

e 2 m^ v 1 ^ 2 ' sm 


7rN\ 


x\ sm 


•TtA^o 


-X2 ■ 


Hence, the obtained Green’s function is correct. 


6. Summary 

In this last section we wish to summarize and point out the main results reached in 
this paper. Our first conclusion is that the separability of the Schrodinger equation of an 
unconstrained system of two particles, coupled or not, generally brakes when boundary 
condition of geometric nature (holonomic constrains) are taken into consideration. The 
main reason of such behavior may be traced back to the ’’coupling” get established 
between the ’’coordinates” as a consequence of the algebraic equations describing the 
same constraints. We have illustrated these point writing the Schrodinger equation of 
two noninteracting particles referred to X c and x, center of mass and relative motion 
coordinate, respectively. By imposing generalized confinement conditions on the system, 
that is limiting the region of motion and assuming impenetrability conditions, one is 
lead to equations relating X c and x which makes, impossible the separation of X c and x 
in the problem. It is remarkable that searching the stationary state of this constrained 
system exhibits the same mathematical formulation one should use to treat the problem 
of a single particle in a plane but confined in a domain whose shape depends on the 
boundary conditions imposed to the two-particle system. Moreover, free particle motion 
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on the plane is studied. The domain of the motion restricted by the impenetrable walls 

forming the square, the rhombus, the triangle and the rectangle (billiards) are considered. 

The billiards are analyzed by studying the time evolution of the covariance of X c and x. 

To this end the Green’s function expressed in terms of the Jacobi $3 function is applied. 
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